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Abstract 

A model  for  the  direct  calculation  of  the  dipole  polarizability  tensors  of  a helix  with  spherical 
core  based  on  the  exact  field  solution  is  introduced.  Motivated  by  this  model,  an  approximate 
one  relying  on  the  dipole  moments  of  the  core  is  developed.  Results  obtained  with  both 
models  are  compared  and  verified  by  measurements. 


1.  Introduction 

Recently,  chiral  materials  consisting  of  complex  inclusions  — for  example  perfectly  electrically 
conducting  (PEC)  thin  wire  helices  with  a dielectric  and/or  magnetic  core  different  from  the  host 
material  — were  proposed  and  investigated  [1].  In  this  paper  the  essentially  exact  model  of  a helix 
with  spherical  core  used  for  the  calculations  in  [1]  is  outlined.  Because  the  governing  equations  are 
lengthy  and  their  derivation  mathematically  involved  only  the  Ansatz  is  displayed  here.  It  reveals 
a problem  concerning  cores  with  both  dielectric  and  magnetic  properties.  The  (approximate) 
solution  of  this  problem  motivates  a second  model  relying  on  a dipole  approximation  of  the  core. 
This  model  is  developed  in  the  following.  Note,  that  cartesian  coordinates  are  used  throughout 
the  whole  paper. 

2.  Model 

Consider  a helix  with  spherical  core  both  centered  at  the  coordinate  origin.  The  current  flowing 
on  the  helix  wire  can  be  calculated  from  a Method  of  Moments  (MoM)  solution  of  the  scattering 
problem  [2],  If  a thin- wire  Galerking  MoM  is  employed,  the  (symmetric)  system  matrix  Z and 
the  exciation  vector  V_  read  for  plane  wave  incidence  (electric  field  E°) 

Zmn  = jufi  f f t(u)(Gi  + g?)t' (u')Bm(u')Bn(u)  du' du 
JuJu' 

Vn  = f f(tt)(/  + £(u))E°  exp{— jfcr}  Bn(u)  du  . 

Ju 

Here,  J is  the  unity  tensor,  G°  stands  for  the  free  space  dyadic  Greens  function,  G*  describes 
the  scattering  correction  caused  by  the  spherical  core  [3],  t (£')  is  the  tangential  vector  on  the 
wire  in  the  observation  (source)  point,  and  B stands  for  the  basis/testing  functions.  The  term 
that  describes  the  plane  wave  scattering  of  the  spherical  core  can  be  derived  from  the  scattering 
correction  G“  by  plane  wave  normalization  as  described  in  [4,  5];  it  is,  in  essence,  identical  to  the 
Mie  solution  for  the  scattering  of  a plane  wave  by  a sphere  [6,  7]. 

The  total  field  scattered  by  this  helix  with  core  can  be  split  into  four  parts,  each  corresponding 
to  an  individual  multipole  series.  These  parts  are:  1)  The  field  scattered  by  the  wire  caused  by 
the  incident  plane  wave,  2)  the  field  scattered  by  the  core  caused  by  the  incident  plane  wave,  3) 
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the  field  scattered  by  the  wire  caused  by  the  diffraction  field  of  the  core,  and  4)  the  field  scattered 
by  the  core  caused  by  the  diffraction  field  of  the  wire. 

Following  the  procedure  introduced  in  [8],  the  incident  plane  wave  in  (1)  can  be  expanded  into 
a Taylor  series  and  the  different  terms  arising  can  be  separated  into  their  electric  and  magnetic 
origin.  This  procedure  gives  results  identical  to  those  obtained  with  the  method  of  counter- 
propagating  waves  [9],  but  removes  the  information  on  the  direction  of  incidence  analytically 
instead  of  using  different  angles  of  incidence  to  separate  the  polarizabilities.  The  Taylor  series 
expansion  produces  unambiguous  results  if  the  spherical  core  is  either  purely  dielectric  or  purely 
magnetic.  If  it  possesses  both  properties,  the  tight  coupling  of  magnetic  and  electric  fields  within 
the  core  prohibits  an  exact  solution.  One  way  to  circumvent  this  problem  is  the  following:  The 
inclusion  under  investigation  is  intended  to  form  the  basic  building  block  of  a chiral  material; 
thus,  the  inclusion  itself  and  all  of  its  constituents  must  be  electrically  small.  Then,  electric  and 
magnetic  effects  are  decoupled  nearly  entirely  (static  limit)  and  can  be  treated  separately  by 
applying  the  series  expansion  of  the  excitation  vector  V_  (see  above). 

This  idea  leads  directly  to  a further  simplification  of  the  problem:  As  the  core  is  electrically 
small,  only  the  first  term  of  the  infinite  series  defining  Ga  contributes  significantly  [3].  This  term 
represents  the  dipole  contributions  to  the  scattered  field.  This  means  that  for  the  description  of 
the  interaction  between  wire  and  core  the  latter  can  effectively  be  replaced  by  its  electric  and 
magnetic  dipole  moments  in  the  origin.  The  expression  for  the  total  electric  field  together  with 
the  boundary  condition  on  the  PEC  wire  surface  then  reads 

Zmn  =ju  [ [ (/x o t_G°  i - uPfilxgPi 0,  u')t'-YV  x £°(0,  u')  tj)  Bm{u)  Bn(u')  du'  du 
Ju  Ju' 

Vn  = 

A = tg°(u,0)^e  , Y = tVxg°(u,0)gm  , 

where  gc  and  gm  are  the  symmetric  electric  and  magnetic  dipole  polarizability  tensors  of  the 
core,  respectively.  To  ensure  reciprocity,  the  identities  t(u)  = t'(u)  and  r(u)  = rf(u)  must  be 
fulfilled  when  considering  Vn  as  derived  in  [8]. 

Together  with  the  abbreviations 

r = 2 zmn  Bn(u')  , (/(u,  u')>  = J J f(u,  V.')  dll'  du 


J (tE^exp{-jkr} 


T -/V  m 


(2) 


the  expansion  of  the  plane  wave  contained  in  (2)  now  leads  to  the  following  dipole  polarizabilities 
of  the  inclusion  (the  numbering  refers  to  the  different  scattered  fields  discussed  above): 

Stuart'),  2e1m  = -?ar(r'xt')) 

JUJ  L 

a™e  = ^<(r  xt)  r t')  , g™71  = - ^((r  x t)  r (r'  x t')) 


=2  = = » £*2 m = £m 


ge3e  = -jutioitrx!) , 


ge3m  = -{trv') 


,2„2 


g™  = —^—((r  x t)  rx!)  , = -^((r  x t)TY!) 


=3 


2 


2 
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g - = -jum&rt) - iJuKxrr) , g«m  = -^(iry  x ())  - ^Vo(irr') 

g7  = (yrt-) + w2w(y  rr> , gT”  = -^<yr  (/  x t')>  - jw^ry') 

Note,  that 

X = iG°(«,0)^e  = ^eg°(0,«)t , Y = tV  x^(u,0)f  = amVx^(0,u)i  . 

Summing  up  all  these  individual  contributions,  the  polarizability  tensors  of  the  whole  inclu- 
sion are  obtained.  It  can  be  proved,  that  the  resulting  tensors  represent  a reciprocal  inclusion, 
i.e.  aee  = (aee)T,  ame  = -(aem)T,  amm  = (amm)T- 

3.  Verification 

If  the  core  is  purely  magnetic  or  purely  dielectric  the  approximate  model  can  be  directly  compared 
to  the  exact  solution.  Furthermore,  results  derived  from  the  system  matrix  Z_  of  eq.(l)  — for 
example  the  resonance  frequency  of  the  inclusion  — are  exact,  even  for  cores  with  both  dielectric 
and  magnetic  properties.  Thus,  these  quantities  can  be  compared  to  the  approximate  ones 
obtained  from  eq.(2)  for  a partial  check  of  the  accuracy.  All  tests  performed  showed,  that  the 
results  calculated  from  the  exact  and  the  approximate  model  are  in  very  good  agreement.  As  an 
example,  the  change  in  resonance  frequency  due  to  a spherical  core  (diameter:  2a  = 3 mm,  helix 
dimensions:  diameter  4 mm,  height  4 mm,  three  turns)  with  magnetic  and/or  electric  properties 
is  displayed  in  Fig.  1.  As  can  be  seen,  the  errors  occurring  in  the  curves  for  either  purely 
dielectric  or  purely  magnetic  cores  add  up  for  cores  with  both  properties.  As  expected,  the  error 
increases  with  electrically  larger  cores.  As  a consequence  of  the  neglected  electric  quadrupole 
moment  this  is  more  pronounced  for  dielectric  cores.  Although  for  er,  [ir  > 5 the  cores  are  no 
longer  electrically  small  (0.2  < ka  < 0.91),  the  results  are  in  good  agreement,  too  (maximum 
difference:  0.04  GHz). 

Fig.  2 displays  the  scattering  parameter  Sn  of  a helix  (dimensions  as  above)  with  a spherical 
dielectric  core  ( eT  = 4.7  — j 1.6,  diameter:  3 mm).  It  was  measured  using  the  method  described 
in  [10]  and  calculated  using  the  approximate  model.  Apart  from  a small  shift  in  the  resonance 
frequency  (see  [11]  for  explanations)  the  results  fit  quite  well,  in  particular  if  one  takes  into 
account  the  measurement  accuracy  [10]. 

4.  Conclusion 

The  approximate  model  for  helices  with  spherical  core  developed  above  produces  meaningful 
results,  both  theoretically  (reciprocity)  and  numerically.  Although  the  model  has  its  origin  in 
the  exact  solution  of  the  helix/spherical  core  problem,  it  contains  no  restrictions  regarding  the 
geometry  of  the  core  or  the  wire.  The  inclusion  should  only  be  electrically  small  so  that  higher 
order  multipole  moments  of  the  core  can  be  neglected.  Thus,  the  model  should  be  able  to  handle 
other  than  spherical  cores,  provided  their  electric  and  magnetic  polarizability  tensors  are  known, 
either  analytically  or  from  numerical  calculations. 
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Figure  1:  Change  in  resonance  frequency  due  to  a spherical  dielectric 
and/or  magnetic  core.  ( ) Approximate  model,  ( ) exact  model. 
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Figure  2:  Scattering  parameter  Su  of  a helix  with  dielectric  core.  Mear 
sured  ( ) and  theoretical  ( ) response  (approximate  model). 


